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' Abstract 

It is argued that the gluon density introduced by A.Mueller [5] in relation to the 
^vq ' saturation phenomenon is different from the standardly defined one, which turned out 

, i , to have a soliton-like form in the numerical calculations. It is shown that A.Mueller's 

densities obtained from the latter possess the same saturation properties as originally 
, predicted. 

m : 

■ 1 Introduction 

> : 

The idea of parton saturation in large nuclei has recently obtained much attention [1-5]. 
In particular A.Mueller has discussed in [5] the behaviour of the quark and gluon densities 
for different values of momenta in a large nucleus. In his treatment the interaction of the 
probing current and the nucleons of the target was taken in the Glauber form, with the 
parton-nucleon amplitude taken in the two-gluon exchange approximation. The evolution 
of the gluonic distribution in x has not been studied. The result of these study is that the 
behaviour of both the quark and gluonic distribution in the momentum I is determined by a 
ft, ■ cerain critical "saturation momentum" defined as 

J§ ■ n 2 8ir 2 a 8 N c 

• • . Qs = C N 2_l AT \ h )i C 1 ) 

• i-H , 

^ ■ with T(b) the standard nuclear profile function and c essentially a constant. At high momenta 

^ . I 2 >> Q 2 the distributions retain their perturbative form 
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d(xq(x,Q 2 )) = N c Q 
d 2 bd 2 l ~ 6vr 4 I 2 

and 



(2) 



d(*G(x,Q 2 )) _ Nl^lQl 1 2 2 

d 2 bd 2 l ~ 47T* I 2 x tOT 1 >>Qs ' {6) 
for the quark and gluon distributions respectively. However as the momentum I diminishes 
both distributions do not infinitely grow (as would follow from (2) and (3)) but saturate at 
a certain finite value. In fact at I 2 << Q 2 A.Mueller found 

d(xq(x,Q 2 )) = jVe_ 

d 2 bd 2 £ 2vr 4 1 ' 

d(xG{x,Q 2 )) j^l 1 {oTi 2 <<0 2 ( 5 ) 

d 2 bdH 47T 4 lU X 1 <<Qg - 

Independently of the saturation problem, in the last year the interaction of a probe with 
a large nucleus at small x was studied in the perturbative QCD with reggeized gluons. It 



and 
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reduces to summing fan diagrams constructed of BFKL pomerons with a splitting triple 
Pomeron vertex [6-9] . In [9] the resulting equation was solved numerically, which enabled us 
to calculate the corresponding gluon distribution. Its features turned out to be somewhat 
different from what has been found by A.Mueller. Although at large I the gluonic distribution 
indeed falls down similarly to the perturbative result and to what has been found in [5], it 
also falls down to zero at small I 2 . The form of the found gluonic distribution as a function 
of log I is close to a Gaussian with a constant height and a center moving with a constant 
velocity to higher I as the rapidity y increases . Thus in the (log l,y) space the gluonic 
density turns out to be a soliton wave. 

This unexpected behaviour of the gluonic density caused some mistrust due to the seem- 
ingly sound base of the analytic derivation of A.Mueller and a purely numerical character of 
the calculations in [9]. Here we intend to study the reasons for the difference of these results. 

One might think that this difference is due to a complicated gluonic interaction introduced 
in the approach pursued in our numerical calculation and absent in the picture of A.Mueller. 
However the answer seems to be much simpler: the gluonic densities studied in both papers 
are in fact different quantities. While the gluonic density numerically calculated in [9] is 
more or less standard, the one introduced by A.Mueller is a new quantity. We show that 
for this latter density our numerical results indeed lead to the same saturation as in [5]. 
Moreover, A.Mueller's quark density found from our solitonic gluon density also exhibits the 
same saturation properties as found in [5]. 

Besides showing that there is no contradiction between [5] and [9] at all, this confirms 
that the gluonic interaction at small x and the resulting very sophisticated dynamical picture 
do not actually influence the saturation properties of A.Mueller's densities, as conjectured in 
[5]. 

Of course there remains a question, which is the "correct" parton density. As we argue 
in the last section, we beleive that the densities introduced by A.Mueller do not in general 
describe partonic distributions in the target. 



2 Gluonic distributions: standard and of A.Mueller 

We begin with the definition of the gluonic density adopted in our numerical calculations, 
which we call "standard". It can be taken from the standard expression for the structure 
function of the nucleus in its terms (see e.g [10 ]) 

F 2 (x,Q 2 ) = J d 2 bd 2 r[ PT (r)+p L (r)]f(x,r,b), (6) 

where pt,l are the transverse and longitudinal colour densities created by the incoming photon 
and / is expressed via the double density of gluons in momentum and impact parameter: 

f(x rb)= f d " k d{xG{x,k 2 ,b)) / _ fer x / _ fer x 

In [9] it was found that, up to a factor, fuction / coincides with a sum of BFKL fan diagrams 

f(*,r,b) = J^-Jj>{x, r ,b). (8) 
Taking a Fourier transform of (7) and neglecting the term proportional to S 2 (k) we obtained 
d(xG(x,k 2 ,b)) N c , N c , ,. 
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where in the coordinate space 

Function </> was found in [9] by numerical solution of the BFKL fan diagram equation. 
Applying k 2 V 2 to it we found the gluon density As mentioned in the Introduction, at any 
given rapidity the density turned out to have the same, roughly Gaussian, shape in variable 
£ = In A;, centered at the point £ = £o(y)> which moves to the right with a nearly constant 
velocity. Approximately the distribution can be described by 

h{k) = h e~ a ^ o(y))2 , (11) 

where ho and a are practically independent of y and Co(u) linearly grows with it: 

ho ~0.3, £ (y) = c + 0.222 mA + 2.23^^y. (12) 

7T 

It can be shown that it is normalized according to 

jH^%,fc,6) = l, (13) 

which condition relates a and /io- The starting position of the center c depends on the initial 
distribution at small rapidities. 

Evidently at a given value of k the density always remains limited, irrespective of the form 
of the initial distribution (and on the atomic number A, in particular). In this sense we have 
saturation as discussed in [1-5]. However with the growth of y the strongly peaked density 
moves away toward higher values of k so that the density at a fixed point tends to zero at 
high values of rapidity. We thus find "supersaturation" : with y — > oo the gluon density at 
an arbitrary finite momentum tends to zero. 

Now we pass to the gluon density introduced by A.Mueller. He related it to a particular 
process, analogous to the 7*-target scattering, in which the electromagnetic current is changed 
to a "gluonic" current 

j(x) = —F^F^. (14) 

This current generates a colourless pair of gluons, which according to [5] interact with the 
target very much in the same way as a qq pair generated by the electromagnetic current 
(see Fig. 1, diagrams b and c). The gluon density is then defined as a contribution to the 
cross-section from a given value of the transverse momentum I of the struck gluon on the left 
part of the diagrams. For the case when the target is just a nucleon, in the limit Q 2 — ► oo 
and in the single gluon exchange exchange approximation, A.Mueller found for his density, 
which we denote xGm(x,Q 2 ,1 2 ) (it generally depends on two different momenta): 

d(xG M (x,Q 2 ,e 2 )) aN c flQ 2 dz d 2 nd 2 r 2 _ u . (n _ r2) ,2(nr 2 ) 2 



d 2 l 2tt 5 Jxfi/Q 



r /w dz grigg -iHr^ ^riY _ 

I 11 V 2 2 / 



(-[ , .-»fc'(ri-rj) _ -ih-n _ ik-r 2 \ ( ^_^_ d(xG(x, k )) . . 

' l+ k 2 dk 2 [ ' 

On the right-hand side there appeared xG(x, k 2 ), which is just the standard gluonic density. 
In the approximations made in deriving (15) it does not depend on x, nor the whole expression 
depends on Q 2 . As a result, the gluon density defined by (15) also depends on only one 
momentum squared £ 2 and superficially looks similar to the standard gluon density. But in 
fact it is a different quantity (factor a s in front of (15) is already sufficient to see it). Its 
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a 



b 



c 



Figure 1: Interaction of a probe with the target. Thick lines show partons (gluons or quarks) 
directly coupled to the probe. Thin lines show exchanged gluons. 



physical meaning has been explained before: it is a probability to find one of the struck 
gluons (the left one in the Fig. 1) with momentum I. Later we shall discuss to what extent 
this quantity probes the gluon distribution in the nucleon. 

We see that once the standard xG is known, the density xGm can be found directly from 
Eq. (15). This exersize is done in the next section. 



3 The soliton density and saturation 

The first thing to do is to generalize xGm to the nucleus target. Having in mind the fan 
diagram picture for the interaction of a probe with the nucleus, on the one hand, and the 
scattering process considered by A.Mueller to define his density, on the other, all we have 
to do is to change the standard density xG on the right- hand side of Eq. (15) to the same 
quantity for the nucleus: 

d(xG M (x,Q 2 ,f,b)) _ aNc dz_ d 2 nd 2 r 2 _ u<ri _ n) 2(nr 2 ) 2 



d 2 bd 2 t 2ir 5 J x e 2 /Q 2 z r\r\ ° ^ r\r\ ^ 

h 4_ p-iHri-rz) _ p-ik-rx _ ik-r 2 \ ( ^_ d(xG) 

1 + ' k 2 d 2 bdk 2 ' [ ' 

Here we suppresed all evident arguments in G on the right-hand side. This will be also often 
done in the following when the arguments are obvious. 

We want to calculate this expression taking for d(xG)/d 2 bdk 2 the solitonic gluon distribu- 
tion in a nucleus obtained by our numerical studies. In principle this can be done rigorously, 
by numerical integration. However, since we are interested in the qualitative features of the 
saturation behaviour, we shall adopt here a simpler approach. We shall use the fact that at 
large rapidities the solitonic distribution is strongly peaked at the momena determined by 

k = e & ~ A 2 / 9 exp ( 2 .23^y) (17) 

(which are very high already at reasonable y for not very small a s ). This momentum will 
play the role of Mueller's saturation momentum Q s . Assuming that k 2 , « Q 2 one can take 
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all the integrand out of the integral over k at point k = k^ except for the density factor. The 
latter, integrated over k 2 gives N c /(2a s -ir 2 ) according to (9) and the normalization condition 
(13). So we are left with the expression 

d(xG M ) _ N 2 1 f JA d 2 nd 2 r 2 m „ !f . (ri „ r2)/ 2(rir 2 ) 2 



l/# d T r2 e-^^^)( 2(r r^ -1) 
x J rir~ rir^ 



d 2 bd 2 £ Air 7 x j , 1 , 2 < ± i 2 

• (1 + e -*(n-r 2 ) _ e -^- r i _ e * fc >^) ( 18 ) 

where one should take k = k^n with a two dimensional unit vector n. 

Calculation of the integrals over r\ and r 2 is simplified by using the identity 

1 (2( n r 2 ) 2 \ 1,__, 2 



1 =7j(ViV 2 r lnrilnr 2 . (19) 



r i r i v r i r i y ^ 

One obtains 

= In i y #[*(/, /) + X(/ + fc, I + A:) - 2X(Z + k, I)] , (20) 

where 

X(/!,/ 2 ) = / d 2 r 1 d 2 r 2 e i ' iri+l ' 2r2 (V 1 V 2 ) 2 lnr 1 lnr 2 = An 2 ^-. (21) 

This leads to 

<*(*Gm) _ iV c 2 , 1 [ , /k 2 -{lk) 2 

Let us consider the two limiting cases studied in [5]. If I » ko then one finds from (22) 

d(xG M ) _ N 2 k 2 1 

d 2 bd 2 i ~ 7T 4 Z 2 X' 1 j 

In the opposite case I 2 «k\ we find 

d(xGu) _ N 2 1 

This behaviour is identical with (3) and (5) (for N c >> 1 and except for factor 4 in (24)). 

So Mueller's density indeed saturates to a finite value at small I 2 irrespective of the 
dynamical picture. The standard density however behaves very differently and for a heavy 
nucleus goes to zero at any fixed value of I as the rapidity gets large enough. 



4 The quark density from the soliton gluonic density 

Following [5] and similarly to the gluon case, one may introduce a quark distribution depend- 
ing on two momenta squared xqM(x, Q 2 ,£ 2 ), as a contribution from the left struck quark to 
the cross-section corresponding to diagrams b and c in Fig. 1. In the two-gluon exchange 
approximation for the interaction, this quantity again results independent of Q 2 at high vir- 
tualities. We postpone discussion of the relation of this limiting distribution in £, xqM(x, £ 2 ) 
to the standard quark density as a function of Q 2 , xq(x,Q 2 ), until the last section. Here 
we shall consider the properties of xqM(x,£ 2 ) which follow from the found standard gluon 
distribution in a large nucleus. 
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The definition of xqM(x,Q 2 ,£ 2 ) is very similar to that of (15). For a single nucleon as a 
target A.Mueller found: 

d{xq(x,Q 2 ,e 2 )) aQ 2 f 1 2 2 / £ 2 (£,£ + k) \ d 2 k dxG 



£ dz[z 2 + (1 - zf] ^ 



d 2 £ 4tt 3 y 1 v /J V£ 2 + e 2 (£ 2 + e 2 )[(£ + k) 2 + e 2 ]J k 2 dk 2 ' 

(25) 
where 

e 2 = Q 2 z(l - z). (26) 

As one observes, the quark density xqu is again defined via the standard gluon density. 

For the nucleus target, as explained before, all we have to do is to take the corresponing 
nuclear density on the right-hand side of (25): 



d{xqu) _ ctQ' 



d 2 bd 2 £ 4vr 3 



Jo 1 ( z >l\p + € 2 (£ 2 + e 2 )[(£ + k) 2 + e 2 ]j k 2 d 2 bdk 2 ' ( n 



At high Q 2 evidently only small z or 1 — z contribute. Then integrating over e 2 ~ Q 2 z 
we find 

d{xg M ) a fd 2 k dxG ( _ (1,1 + k) (I + k) 2 \ 

d 2 bd 2 £ 2tt 3 J k 2 d 2 bdk 2 { (I + k) 2 -I 2 I 2 )' ( ' 

Let us take for xG the soliton-like gluon density found in [9]. As before we take all smooth 
functions out of the integral over k and afterwards integrate over k using (13). We obtain 



d{xq M ) N c f r (1,1 + k) (l + k) 2 



(I + k) 2 - I 2 I 2 

where k = k^n as in the preceding section. 

Again we shall concentrate on the two limiting cases. If I » ko using 

1 ^ (/ + k) 2 _ 1_ ^ 1 2(lk) + k 2 ( 4 (Ik) 2 



(29) 



In - — ~ -7T 1 h 



(I + fc)2 _p p - p \ 2 p 3 ? 4 

and integrating over the angle 4> w e get 

d(xq M ) _ Nc fep 
d 2 fe(i 2 £ ~ 6tt 4 Z 2 ' 



(30) 



In the opposite case I 2 << ko the second term in the bracket of (29) does not contribute and 
we get 

d(xq M ) _ N c 



d 2 bd 2 £ 2vr 4 ' 

This behaviour exactly coincides with A.Mueller results (2) and (4). 



(31) 



5 Discussion 

As we have seen, there are in fact two different gluon densities, which exhibit different be- 
haviour as functions of momenta. The standard one defined via Eq (7) has a soliton-like 
behaviour and vanishes at any given £ as y — > oo. The other, defined in [5], saturates at 
small momenta to a finite value, which turns out to be universal in the sense that it does not 
depend on the complexity of the interaction mechanism. Which of the two densities is the 
"correct" one? 

The definition of the standard gluon density via Eq. (7) seems to us perfectly adequate. 
Its validity is clear from the inspection of the diagrams b and c of the Fig. 1 corresponding 
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to (7). The lower gluonic blob, to which the quark loop is coupled, has to be integrated over 
the "-" component of the 4-momentum k (in a system, in which the target has momentum 
with a large "+" component). Evidently ihe blob then represents just the gluonic density at 
a given "time" x+. 

Interpretation of Mueller's density, in contrast, is not so obvious. First, in the presence 
of the interference diagram c of Fig. 1, one may hardly say that the contribution from the 
left struck gluon to the cross-section probes the gluon distribution in the target. In the 
interference diagram the right struck gluon has a different transverse momentum £ + k, so 
that the interpretation of the whole contribution as a gluon density is doubtful. Then it 
is not clear why the lowest order diagram Fig. la has not been included. Finally, in the 
BFKL kinematics the virtual gluons tend to reggeize. One then expects that the interference 
digram can be neglected and diagram Fig. lb will be absorbed into a to form a BFKL 
pomeron coupled to the target. So in the end one will find only the diagram Fig. la, which 
leads to the standard gluon density as defined through Eq. (7). 

The same doubts refer to the quark density in the definition of A.Mueller. Again, so 
long as one cannot neglect the interference diagram, the interpretation of the contribution 
to the cross-section from the left struck quark with a given momentum as a quark density 
of the target does not seem possible. True, in the DGLAP kinematics, with finite x , the 
interference diagram can indeed be neglected, substituted by the requirement that £ 2 << Q 2 . 
However then one has the transverse momentum ordering and the density found from Eq. 
(25) would be trivial (see [5], Eq. (14)). In the BFKL regime, at small x, the virtual quarks 
also tend to reggeize and the contribution from digram Fig. lc hopefully becomes neglegible 
as compared to a and b, which actually become identical. Then the interpretation of (25) 
as the quark density may aquire sense. However to calculate it one should couple the found 
gluonic distribution to a reggeized quark. We leave this problem for future studies. 
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